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Abstract. This paper is devoted to the study of numerical approximation schemes for 
the heat equation on (0, 1) perturbed by a non-linear rough signal. It is the continuation 
of IE] , where the existence and uniqueness of a solution has been established. The 
approach combines rough paths methods with standard considerations on discretizing 
stochastic PDEs. The results apply to a geometric 2-rough path, which covers the case 
of the multidimensional fractional Brownian motion with Hurst index H > 1/3. 



1. Introduction 

This paper is part of an ongoing project whose general objective is to extend the scope 
of applications of the rough paths method to infinite-dimensional equation, with as a 
target the possibility of a pathwise approach to stochastic PDEs (see [151 El SI ITT])- 
The equation we mean to focus on here is the following: 



yo 



^GL 2 (0,1) , dy t = Ay t dt + J2fi(yt)dxi , t G [0, 1], (1) 



where: 



• A is the Laplacian operator on L 2 (0, 1) with Dirichlet boundary conditions, 

• fi(Vt)(0 ■= fi(Ut(0) for some regular function /■ : R ->■ R, 

• x : [0, 1] — > R m is a geometric rough path of order 1 (see Assumption (Xl) 7 ) or 
2 (see Assumption (X2) 7 ). 

Owing to the results of [5] , we know that the latter hypothesis includes in particular the 
case where a; is a fractional Brownian motion (fBm in the sequel) with Hurst index H > 
1/3. Thus, Equation ((T|) provides in this situation a model that can deal with the long- 
range dependance property at the core of many applications in engineering, biophysics 
or mathematical finance (see for instance [7J [221 l2Sj). It also worth mentionning that 
in the fBm case, the equation can also be handled with Malliavin calculus tools (see 
[2"9"| [2U I2S1 HB]), but for H > 1/2 or for very particular choices of /j only (/; = 1 or 
fi = M). 

The theoretical treatment of ([T]) under its general form has been established in [9] and 
[Sj. More precisely: 

(i) When x is a geometric 1-rough path, it is proved in [9] that ([[)) admits a unique 
global solution for any regular enough initial condition ip, and this is obtained by means 
of an abstract fixed-point argument in a well-chosen class of processes. 
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(ii) When a; is a geometric 2-rough path, the existence and uniqueness of a global solution 
has been shown in [8] via a time- discretization of the equation. 

We will go back to the exact statement of those two results in Sections [3] and HI Let 
us only point out here that in both situations, explicit solutions are rarely known and 
the arguments at the basis of these existence results are not sufficiently constructive to 
provide a representation of the solution. This paper is meant to remedy this problem 
by introducing easily-implementable numerical schemes for the two configurations (i) 
and (ii). The approximation procedure will stem from two successive discretizations, in 
accordance with the classical strategy displayed for Wiener SPDEs (see [IB] or [17]): we 
first turn to a time- discretization of the problem and then perform a space-discretization 
of the algorithm, following the Galerkin projection method. 

The schemes will actually be derived from the theoretical interpretations of ([TJ con- 
tained in [91 IB]. For this reason, let us remind the reader with a few key-points of the 
approach displayed in the latter references: 

• The equation is in fact analyzed in its mild form, namely 

y t = S t il>+ f S t . u dx\ fi{y u ) , t G [0, 1], (2) 
Jo 

where S stands for the semigroup generated by A. This is a classical change of 
viewpoint in the study of (stochastic) PDEs (see [6]), which allows to resort to the 
numerous regularizing properties of S (summed up in Subsection 12.51) . 

• As with rough standard systems, the interpretation of the right-hand-side of (jSJ) relies 
on the expansion of the convolutional integral J* S t - U dx l u fi(y u ), which gives rise to a 
decomposition such as 

t 

St-u dx l u fi(y u ) = Pts + Rtsi (3) 

where P is a "main" term and R a "residual" term of high regularity w.r.t (s,t), 
which is likely to disappear from an infinitesimal point of view. Once endowed with 
this decomposition, the time-discretization is naturally obtained by keeping the main 
term P only between two successive times of the partition: 

yo=^ , y?* +1 = s tk+1 - tk y tk + P tk+ltk , (4) 

with for instance tk = tff = k/M. The reasoning can here be compared with the recent 
approach of Jentzen and Kloeden for the treatment of a Wiener noise ([HI [201 121]): 
in order to deduce efficient approximation schemes, the two authors lean on a Taylor 
expansion of the solution, which indeed fits the pattern given by ([3]). 

• Then, in comparison with the standard case, an additional step has to be performed 
so as to retrieve a practically-implementable algorithm: roughly speaking, it consists 
in projecting the (intermediate) scheme (jlj) onto (increasing) finite-dimensional sub- 
spaces of L 2 (0, 1). We will thus carefully examine how to combine this projection with 
the rough paths machinery. 

To the best of our knowledge, this is the first occurence of (explicit) approximation 
schemes for a PDE involving a fractional noise. The convergence of those schemes will 
hold for any geometric 2-rough path. We hope that the strategy as well as the technical 
arguments displayed in this paper will make possible the approximation of a larger class 



NUMERICAL SCHEMES FOR THE ROUGH HEAT EQUATION 



3 



of rough evolution equations, with for instance a more general operator or a fractional 
distribution-valued noise. For the time being, we cannot handle this task though, just 
because theoretical (global) solutions have not been obtained in those situations yet. 

The article is organized as follows: in Section [21 we first elaborate on the assumptions 
underlying our study. We also introduce the two algorithms that will be brought into 
play and state the main convergence results. Section [3] is devoted to the treatment of the 
above case (i). Only developments of order 1 will be involved in this section, so that the 
scheme can be seen as an adapted version of the usual Euler scheme. In Section HJ we 
will handle the scheme for the situation (ii) , which requires developments of order 2 and 
is thus closer to the well-known Milstein approximation for standard differential systems. 
Finally, Appendix A puts together some technical proofs that have been postponed for 
sake of clarity, while in Appendix B we give an insight into possible implementations of 
the algorithm in the fBm case. 

2. Settings and main results 

2.1. Framework. We focus on the Laplacian operator A on to the Hilbert space B := 
L 2 (0, 1) with Dirichlet boundary conditions. We fix from now on a basis of B made of 
eigenvectors: 

e n(£) := V / 2sin(7rn^) (n £ N*), with associated eigenvalues A n := i^n 2 . 

For any N £ N*, Pjy wm stand for the projection operator onto the finite-dimensional 
subspace Vn ■= Vect {e n , 1 < n < N}. It is a well-known fact that the fractional 
Sobolev spaces are likely to play a prominent role for the study of a stochastic PDE 
(see e.g. [23]): 



Notation 2.1. For any k > 0, we denote by B K the fractional Sobolev space associated 
to (— A) K and characterized by 

oo 

B K = {y £ L 2 (0, 1) : ^A 2K (y") 2 <oo}, (5) 

n=l 

where the (y n ) are the components of y in the basis (e n ). This space is naturally provided 
with the norm 

oo 

\\y\\k = \\(- A ) K y\\ 2 B = J2 x n(y n ) 2 , (6) 

n=l 

and we extend the definition of B K to any k < through the characterization formula 

m- 

2.2. Assumptions. As in [HI [9] , we are interested in the mild formulation of the equa- 
tion, namely 

yt = S t 4>+ [ S tu dxi h{y u ) , t £ [0, 1] , ^ e B, (7) 

where S is the semigroup generated by A and S tu := S t - U - A priori, the equation 
only makes sense for a regular (ie piecewise differentiable) process x. In this context, 
interpreting the rough version of ([7]) means extending the convolutional integral to a 
7-H61der process x, 7 £ (0, 1). For sake of simplicity, we will only consider in this paper 
the case 7 is strictly greater than 1/3, which covers in particular the Brownian motion 
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case. As in the classical rough paths theory, we will also be led to assume, depending 
on the regularity of x, that one of the two following assumptions is satisfied. 

Assumption (Xl) 7 : x : [0, 1] — > ¥L m is a geometric 1-rough path of order 7. In other 
words, x is a 7-Holder process and there exists a sequence of piecewise differentiable 
process (x M ) such that 

u M := M[x - x M - C?([0, 1]; R m )] 0, 
where Af[.; CJ([0, 1]; R m )] is just the usual Holder norm, ie 

\yt - j/ 8 ||r™ 



A%;C7([0,l];ir)] := sup 



0<s<t<l 



t - s 



Assumption (X2) 7 : x : [0, 1] — > W 71 is a geometric 2-rough path of order 7. In other 
words, a; is a 7-Holder process and there exists a sequence of piecewise differentiable 

process (x M ) such that M[x - x M ; CJ([0, 1]; R m )] and the sequence (x 2,M ) of 

Levy areas associated to (x M ), ie 

x^'-" := f <!■>■:! ■ (x™> j - ./•;'-'). /../ 1 8 < t e [0, 1], 

J s 

converges to an element x 2 with respect to the norm 



A%C^([0,1];R^)] := sup 



1 2/ts ||r"».™ 



0<s<t<l t — s 



27 



In brief, 

% := A/> - x M ; C7([0, 1]; R m )\ + A^[x 2 ' M - x 2 ; C 2 27 ([0, 1]; M m ' m )] 0. 



Example: The main process we have in mind in this paper is the (m-dimensional) 
fractional Brownian motion x = B with Hurst index H > 1/3. It has been indeed 
proved in [S] that this process satisfies Assumption (X2) 7 (and accordingly Assumption 
(Xl) 7 ) for any 1/3 < 7 < H, when taking for x M the linear interpolation of x with 
uniform mesh ie 

tk = t " := Ji ' x ^ :=x * fc + M -( t -^)-K +1 -^) ifieMfc+i)- (8) 

Besides, the following sharp estimates have been established in [ID] : for any 1/3 < 7 < 
if, there exists an almost surely finite random variable C 7 such that t> m < C^^log M • 
Mi-H _ fj^jg con trol can then be injected in our final estimates (fl^l) and fTl~5]) so as to 
retrieve explicit (almost sure) rates in this fBm case. Note that Conditions (Xl) 7 or 
(X2) 7 are actually fulfilled by a larger class of Gaussian processes, as reported in [T3] 
or in [T2] . 

As far as the regularity of the vector field / is concerned, it will be governed by the 
additional condition (k is a parameter in N): 

Assumption (F)k'- for every i G {1, . . . ,m}, fi belongs to the space C fc,b (M;M) of 
/c-time differentiable functions, bounded, with bounded derivatives. 
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2.3. Schemes. In order to introduce the two schemes we intend to study, let us define, 
for any piecewise differentiable process x : [0, 1] — > M m , the following operator-valued 
processes: for every i, j = 1, . . . , m, for any s < t G [0,1], 

Xff := f S tu dx\ , X™> ij : f S tu dx\ {.r'„ - xj), (9) 

J s J s 

We suppose in addition that either Assumption (Xl) 7 or Assumption (X2) 7 is satisfied, 
for some parameter 7 G (0, 1) and some regularizing sequence {x M ), and that Assumption 
(F)! holds true, so that f- is well-defined. With those conditions in mind, here is the 
two schemes that will come into play in the sequel: 

Euler scheme: Dq I,N = Pn^ and 

M.N n M.N , v x M A D tl M,N\ /-, n \ 

Vt^ = S tk+ltk y tk ' +X tk+ ' itk P N f i (y tk ' ), (10) 

where t k = tf = ±. 

Milstein scheme: Dq I)N = Pn^ and 

M.N a M.N . v x 2M ,i j-, t ( M,N\ . v x 2M x 2M ,ij n ( t i , M,N\ / n t / M,N\\ 

Vt k U = S t k+i t k y tk + X tk+1 i k P N fi{y tk ) + X tk+itk J P N [f^ ) ■ {PNfAVtu )) ' 

(11) 

where tk —t^f — -Mj and the notation (j) ■ ip stands for the pointwise product of functions, 
ie (0-^(0 := 0(0^(0- 

Remark 2.2. The name we have given to the schemes is just a reference to the classical 
algorithms for standard stochastic differential equations. It indicates that ( fTu]) involves 
developments of order one only, while (llip appeal to second-order terms. 

Remark 2.3. When x M is the linear interpolation of x given by (jHJ), the two sequences 
of operators Af 'I , Af ' tJ that intervene in the schemes reduce to 

M ■ ■ ■ f tk + 1 

K + k=M-{x\ k+i -x\ k ). / S tk+lU du (12) 

Jt k 

<rC' ij = M 2 • (xl k+i -x\ k ). (xl k+i - 4) • T +1 ^ fe+lU rf M (it - t fc ). (13) 

Consequently, in this case, the computations of Formulas (|T0i) and f[TT]l only require the a 
priori knowledge of the successive increments Xt k+1 —Xt k , which makes the implementation 
of the algorithms very easy, as we shall see in Appendix B for the fBm case. 

2.4. Main results. We are now in position to state our main convergence results. As 
in the standard rough paths theory results, we make a clear distinction between the case 
7 > \ and the case 7 G (|,|). Theorem 12.41 (resp. Theorem 12.51) will be proved in 
Section [3] (resp. Section HJ) . 

Theorem 2.4. Let 7 G (|, 1) and suppose that Assumptions (Xl)^ and (F)i are sat- 
isfied. Let also 7' G (max(l — 7, |) and ^ G -By. JTien t/iere exists a function 
C : (M + ) 2 — > M + bounded on bounded sets such that if y is the solution of ffl) with initial 
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condition ip ( see Theorem \3. 5\) and y M > M is the process generated by the Euler scheme 
(ED with M = N, 

keo,..,M ^ 
where we have used the shortcut \\x\\^ := J\f[x;Cj([0, l];R m )]. 

Theorem 2.5. Let 7 £ (|, |) and suppose that Assumptions (X2)^ and (F)% are sat- 
isfied. Let also 7' £ (1 — 7, 27] and ip £ Hy. Taen /or every parameters 

< /3 < inf ^7 + 7' - 1,7 -7' + , 0<A<7 + 7'-l, 

there exists a function C = Cp t \ : (M + ) 2 — > R + bounded on bounded sets such that if y 
is the solution of |7|] in By with initial condition ip (see Theorem \4-5\ ) and y M > N is the 
process generated by the Milstein scheme 



sup 

ke{o,...,2M} 



vtif-yf/K < c (||^|| By , ||x|| 7 ) 1 1| v - PniPWb,, + v 2 m + + 

where we have used the shortcut ||x|| 7 := Af[x; Cj{[0, 1]; R m )] + A/"[x 2 ; C 2 27 ([0, 1]; 1 



1 1 



Remark 2.6. The particular choice N = M in Theorem 12.41 has only been made so as to 
get a nice expression for the final estimate ffT4l . Nevertheless, it is not hard to obtain a 
more general result with possibly different N, M, following the arguments of Section |3j 

Remark 2.7. As we shall see in Section HI the use of dyadic intervals in the Milstein 
scheme (TTTj) is justified by the need of a decreasing sequence of partitions in the patching 
argument of Proposition 14.91 However, our convergence result can probably be extended 
to any sequence of partitions whose meshes tend to 0, at the price of more intricate local 
considerations in the proof of the latter proposition. 

2.5. Tools of algebraic integration. Before going further, let us draw up a list of the 
properties at our disposal as far as the fractional Sobolev spaces B K and the semigroup 
are concerned (the proof of those classical results can be found in [2], [25] or |27j): 

• Sobolev inclusions: If k > 1/4, B K is a Banach algebra continuously included in the 
space L°°([0, 1]) of bounded functions on [0, 1]. 

• Projection: For all < k < a and for any if £ B a , 

\\tp - Pnv\\b k < A^ (a ~ K) |M| BQ . (16) 

• Contraction: For any n > 0, S is a contraction operator on B K . 

• Regularization: For any t > and for all — 00 < k < a < 00, St sends B K into B a and 

\\SMK <c Q , K t- (Q - K) |Mk. (17) 

• Holder regularity: For all t > 0, a > and for any ip £ B a , 

\\S t ip - <p\\ B < c a t a \\<p\\ Ba , \\AS t ip\\ B < c a r 1+a y\\ Ba . (18) 
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• Composition: if k G [0,1/2) and / G C 1,b (IR; M) (see Assumption (F)^ for the latter 
notation), 

\\mh K <c K , f {i + y\\ BK }, (19) 

while if k G [1/2, 1) and / G C 2 - b (M;M), 

\\m\\B K <c Kj {l + Ml K }, (20) 
where, in both cases, f((p) is understood in the sense of composition, ie /(</?)(£) : = 

/(</?(£))■ 

• Pointwise product: if k G [0, 1/2) and ^,-0 G i3 K n L°°([0, 1]), 

UP "01k. < ^{IMU^IMk + IMklHU°°}> (21) 
while if <^ G £>_ K , ip G -B^, with k > and a > max(K, |), 

h-^\\B^<c K>a \\<f\\ B _ K U\\ Ba . (22) 
Remember that (p ■ ip is understood as a pointwise product, ie (ip ■ V0(O — ^(OVKO- 

With these properties in hand, the rough paths treatment of Equation (126]) is based on 
the controlled expansion of the convolutional integral 

1 S tu dx l u f t (y u ). (23) 



In order to express this expansion with the highest accuracy, we provide ourselves with a 
few tools ans notations inspired by the algebraic integration theory for standard systems 
(see P3|). 

Notations. For k G {1, 2, 3} and for any interval / C [0, 1], denote 

S fc (I) := {(ti, -..,**) G I k : ti > ... > t k }. 
Then for all processes y : I —> B and z : S 2 (I) —> B, we set, for s < u < t G /, 

(Sy)ts ■= Vt-Vs , (Sy)te ■= (Sy)ts - a ts y s , (24) 

{fifties '■= %ts ~ z tu ~ Stu z us, (25) 

where a ts := S ts — Id. 

To give an idea on how those operators arise from the handling of (JZJ), let us observe for 
instance that the variations of the solution y are governed by the equation 

rt ps r-t 

(tiy)ts = / S tu dx l u fi(y u ) + a ts \ S su dx l u fi(y u ) = / S tu dx l u fi(y u ) + a ts y s , 

Js J0 Js 

and can thus be equivalently written as 



Vo = ^ , (Syh = J Studxlfiiyu). (26) 

Let us also observe, in this convolutional context, the following elementary properties, 
that we label for further use: 

Proposition 2.8. Let y : [0, 1] — Y B, z : S% — > B, and let x : [0, 1] — > R a regular 
process. Then it holds: 

• Telescopic sum: 5(5y) tus = and (8y) ts = YaZq Stt i+1 {o~y)t i+1 ti for any partition 
{s = t < ti < . . . < t n = t} of an interval [s, t] of [0, 1]. 



8 



AURELIEN DEYA 



• Chasles relation: if J ts := J* S tu dx u y u , then 5 J = 0. 

• Cohomology: if 5z = 0, there exists h : [0, 1] — > B such that 5h = z. 

On top of those algebraic considerations, if one wants to measure the regularity of the 
terms involved in the expansion of f S tu dx l u fi(y u ), one is led to introduce the following 
suitable semi-norms, that can be seen as generalizations of the classical Holder norm: 
thus, if y : [0, 1] — > V, z : S 2 — > V and h : S 3 — > V, where V is any Banach space, we 
will denote, for any A > 0, 

Af[y;C?([a,b];V)]:= sup MM^X , ^^([a, b]; V)} := sup \\y t \\ v , (27) 

a<s<t<b \t — S\ te[a,b] 

M[z-C^b}-V)}:= sup , M[h-C${[aM;V)}:= sup 1^.(28) 

a<s<t<b \t — S\ a<s<u<t<b \t — S\ 

Then C^([a, b];V) naturally stands for the set of processes y : [0, 1] — > V such that 
Af[y;C£{[a,b];V)] < oo, and we define C?([a, 6]; V), C 2 A ([a, b]; V) and C^([a,b];V) along 
the same line. With these notations, observe for instance that if y e C% ([a, 6]; £(V, W)) 
and z e C2([a,b];V), the process h defined as h tus := yt u z us (s < u < t) belongs to 
C^([a,b];W). ' 

When [a, b] = [0, 1], we will more simply write C£(V) := C£([a, b]; V). 

The following notational convention also turns out to be useful as far as products of 
processes are concerned: 

Notation 2.9. If g : S n — > C(V, W) and h : S m — > W , then the product gh : «S n+m _i — > 
W is defined by the formula 

(gh)t 1 ...t m+n - 1 '■= 9t 1 ...t„ht n ...t n+m _ 1 - 

With this convention, it is readily checked that if g : <S 2 — > C(B K , B a ) and h : S n — >■ B K , 
then S(gh) : S n+ \ — > B a is given by 

5(gh) = (5g)h-g(5h). (29) 



To end up with this toolbox, let us report one of the cornerstone results of [J5], which 
will allow us, in Section HI to cope with the high-order terms poping out of the expansion 
of 



Theorem 2.10. Fix an interval I C [0, 1], a parameter k > and let fi > 1. For any 

h G C3 (/; B K ) H Im5, there exists a unique element 

Ahen aeM C^ a (I;B K+a ) 
such that S(Ah) = h. Moreover, Ah satisfies the following contraction property: for all 

a e [0,11), 

A/lA/i;Cr Q (/;i3 K+a )] < c aifl M[h;C^I;B K )}. (30) 
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3. Young case 

This section is devoted to the proof of Theorem 12.41 Consequently, we fix from now on 
the two parameters 7 £ (5, 1) and 7' £ (max(|, 1 —7), |), as well as the initial condition 
ip £ By. Under Assumptions (Xl) 7 and (F) 2 , the convolution integral ( l2"B"j) can be 
extended to x via a first-order expansion. To do so, observe that if x is a piecewise 
different iable process, one has, for any £>-valued differentiable process z, 

1 S tu dxl z u = Xffz s + A ts (X*' l 5z), (31) 

where X x ' 1 is the operator- valued process defined by Indeed, if we denote 

Jts ■= / S tu dx l u z u - X*fz s = / S tu dx u (5z)us £ Cf (£), 

J s J s 

one has, with the help of Theorem 12. 10[ 5(J — A(X x,l 5z)) = 0, hence, owing to Proposi- 
tion J — A(X x ' l Sz) — 5h £ C\{B\ which easily entails Sh = (use the telescopic-sum 
property of Proposition ^. 8p . One can then rely on the following natural extension result: 

Lemma 3.1 ([!]). Under Assumption (Xl) 1; the sequence of operator-valued processes 

Stu dx u ' l 

converges to an element X x,t with respect to the topology of the spaces C2~ X (C(B K , B K+ \)) 
(X £ [0,7),k £ R) and Af[X x ^;C^~ x {C{B K ,B K+x ))} < c« )A ||x|| 7 , as well as 

M[XXti _ xxM>i _ C 7-a (£(b ^ Bk+x))] < Ck xUm (32) 

Moreover, commutes with the projection and satisfies the algebraic relation 

5X X > 1 = 0. 

Remark 3.2. The underlying topology of this convergence result is of course closely 
related to the properties of the semigroup recalled in Subsection 12.51 In other words, 
the fact that X x,t £ C2~ X (£(B ti , B K+ \)) for k £ R, A £ [0,7), is a consequence of the 
regularizing property (|T7|) . 

Remark 3.3. Through the continuity result (I32j) . one can see that the process X x only 
depends on x and not on the particular approximating sequence x M . This comment also 
holds for the forthcoming Lemma 14.11 

Once endowed with X x , it is readily checked that the right-hand-side of ( |3T|) can also 
be extended to a less regular process z, which provides the expected interpretation: 

Proposition 3.4 ([9]). Under Assumption (Xl) lt we define, for any process z = 
(z 1 ,..., z m ) such that z i £ C^{B K ) n Cf (B) with 7 + k > 1, the integral 

J ts {dx z) := X x t tz\ + A ts (X^Sz 1 ) . (33) 

Then: 

• J{dxz) is well-defined via Theorem \2.1(A It coincides with the Lebesgue integral 
J s S tu dx u z l u when x is a piecewise differentiable process. 

• The following estimate holds true: 

Af[J(dxz);Cj(B K )} < c\\x\\t{M[z;C?(B k )) +N[z;C?(B)}} . (34) 
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3.1. Previous results. The main result of [9] for the Young case is summed up by the 
following statement: 

Theorem 3.5 ([9]). Under Assumptions (Xl)^ and (F) 2 , Equation ^jj, interpreted 
thanks to the previous proposition, admits a unique solution y in CJ {By), and the fol- 
lowing estimates hold true: 

M[y-A^)] + M[y;C?(By)\ < C (|M| By , |M| 7 ) , (35) 

for some function C : (IR + ) 2 — > IR + bounded on bounded sets. Morever, if y (resp. y) is 
the solution of |7|) associated to a process x (resp. x) that satisfies Assumption (Xl)^, 
with initial condition ip (resp. ip) in By, 

M[y - y- C?(By)] < C X -^{U ~ + \\x- 5|| 7 } , (36) 

with c x,xrt,j '■= C'CII^IIt' II^IIy> IMIb v > ll^lby); f or some function C : (IR+) 4 M + 
bounded on bounded sets. 

Remark 3.6. It is worth noticing that (1351) and (134)) implies in particular 

M[y;Cj(By)} <c^ x . 
Indeed, since y is solution to the system, one has 

\\($y)tshy < \\Jts(dxf(y))\\By 

< c x \t-sr {M[f(y); Cl{By)} + JV[f(y); C?(B)}} . 

Then, thanks to and (EES}, it holds Af[f{y);C%{By)} < c{l + N[y;C%{By)}} and 

JV[f(y);C?(B)} < cM[y;Cl'(B)} < c {M&C^By)} + M[y;Ct (By)}} . 

The continuity result (13"6"j) provides us with a control over the discretization of the 
driving signal x. This is the first step towards Theorem 12.41 

Notation 3.7. For any M £ N*, we denote by y M the Wong-Zakai approximation 
associated to x M (with the same initial condition ip), or otherwise stated the solution to 
Equation |7|) when x is replaced with its interpolation x M . 

Corollary 3.8. With the above notations, there exists a function C : (IR + ) 2 — > M + 
bounded on bounded sets such that, for any M £ N* ; 

sup \\y tk -yfihy < C(\\x\\^, \\i>\\By)u M - (37) 

fc€{0,...,A/} 

3.2. A uniform control. The second step of our reasoning consists in controlling the 
process y M ' N generated by ( fTUj) . uniformly in M and N. To do so, let us first extend 
yM,N on ^ through the formula: if t £ tfc+i), 

M,N M,N . v x M ,i j-, r i M,N\ / 00 \ 

Vt ■= S Uk y tk + x u k p Nfi{y tk )• (38) 

Now observe that by setting r^' N := A ts i^X^' ,l 5PNfi(y M ' N )^ , one can write, for any 
k £ {0,...,M- 1}, 

M,N a M,N . i %M r> r / M,N\ M,N / Qn \ 

y ik+1 = S tk+1 t k y tk + S tk+lU dx^ P N fi{y u ' )-r tk ' +ltk . (39) 



NUMERICAL SCHEMES FOR THE ROUGH HEAT EQUATION 



11 



Extending the expression to all times s < t gives rise to the two formulas: 
Lemma 3.9. If t p < s < t p+ i < . . . < t q < t < t q+ \, then 



with 



(Sy M ' N h 



M,N\ _ ,M,NA 
iJts 1 



M,N q M,N , \ ^ Q M,A 

Vts ■— r tt q ~ dt s r stp + 2^ * ifc +i r *fe+i 



q-1 

£ 

k=p 



N 



while if t p < s < t < t p+ \, 



(Sy M >% s = X?:''P N f^ N ). 



(40) 



(41) 



(42) 



Proof. Formula (142 p is a straightforward consequence of the relation SX xM,i = 0. As for 
ID]) , it follows from the association of ( 139]) and the telescopic-sum property contained 



in Proposition 12. 8[ which gives here 



(Sy M ' N h 



s ttk+1 (ty M > N )t k+1 t k + (h M ' N )u q + S ttp+1 (5y M ' N )t p+ls 

k=p+l 

^ S tu d X - M P Nfl (y^)- £ S ttk+1 r%f lth 

P+ 1 k=p+l 

* S tu dx^ M P N f l (yi I ' N )-r 

1 



+ 



tt a 



+ s ttp+1 (Sy M ' N ) tp+ls (43) 



Since 



(k M > N h +ls = (k M ' N )t P+ltp - S tp+ls (k M > N U P 



t p +i 



S tp+lU dx^ M P N f t (y* 



,M,N\ | M,N 



-5' 



t P +is 



r S su dx^ P N fi(y. 

J tri 



M,N\ _ M,N 
I ' st v 



it suffices to inject (jUj) in to get QMS). 



(44) 



□ 



We are going to lean on the two expressions fT4T3]) and (j42]) in order to establish the 
expected uniform estimate: 

Proposition 3.10. There exists a function C : (R + ) 2 — > R + bounded on bounded sets 
such that for every M, N e N, 



Af[y M > N ; C°([0, 1], By )] + Af[y M > N ; #([0, 1], By)} <C(\\x\ 
where y M ' N is extended on [0, 1] through Formula / TF5]) . 



7' \\W\\Byh 



(45) 
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Proof. For sake of conciseness, denote here 

U[y M,N. d 0,Y (/)] ;= U[y M,N. C (J) By)] + ^M,*. C V (J) By)] . 

With this notation in hand, we will actually prove the following assertion: there exists 
a time T = T (||x|| 7 ) > and a sequence of radii Ri = i?/(||x|| 7 , ||^||b /) such that for 
any /, 

Af[y^ N ;C^'([0,lT })]<R, 
For I = 0, take Ro := ||^||b /• Now assume that the property holds true for I, and let 

s,t e [o,(/ + i)T ]. 

1 st case: s,t G [IT , (/ + 1)T ]. 

I s * subcase: t p < s < t p+ i < . . . < t q < t < t q+ \, with \t — s\ > -jj. Then, from (TJO]), 



J s 

Owing to the estimate (IMI) (applied to x = x M ), one easily deduces 

S tu dx^ P N f i {y^ N )\\ Bi> <c x \t- sfT^' {l +Af[y M > N ;C°S([0, (I + 1)T ])]} 

Besides, thanks to the contraction property (1301) of A, one gets 

ll^lls <c x \t- s\^' {l + Af^; C^'([0, (/ + 1)T ])]} , 



as well as 
Thus, 



K^lk, < c, |t - ,r {l+J^VflO, (J + 1)T ])]} . 



9-2 

II -A4",JV,tt|i ^ || MJVii . n AfJVii , n A/, AT n U j- |-Y II A/, AT n 

WVts Hb v < Hr^ ||b v + ||r stp ||s v + ||r M ;_ 1 ||B T , + 2^ I* - <fc+i| IK +1 t fe b 



< c4i+Ar[/^;C^'([o,(/ + i)r ])]). 



9-2 



< c, {l +Af[y M ' N ;C°S'([0, (I + 1)T ])]} j|t - s| 7 + '^Ij 

< c x |t - «r {l+M[y M ' N ;C^'([0, (/ + 1)T ])]} . 

r d subcase: t p < s < t < t p+1 . Then (8y M,N )ts = X? s M,i P N fi(y™' N )i so that 
||(<fy M,A %lk, < c, |t - S | 7 (l +Af[y M > N ;C°S\[0, (I + 1)T ])] j . 



3 rd subcase: t p < s < t p+ \ < t < t p+ 2 with \t — s\ < 1/M. Just notice that 

||(fy M '*%||s y < \\(k M ' N )u p+1 h yl + \\(k M ' N )t P+1 s\\B^„ so that we can go back to the 
second subcase. 
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Conclusion of the 1 st case: 

M[y M ' N ;CU[lT , (I + 1)T ])] < c x T^ {l + Af[y M ' N ; C^'([0, (I + 1)T„])] j . 



2 nd case: s < IT Q < t < (I + 1)T . One has \\(Sy M ' N )ts\\B yl < II (k M ' N )t,m || By + 
\\(Sy M ' N )iT , s \\b an d s °) owing to the recurrence assumption, 

\\(6y M > N ) ts \\ By , < \t-s\~*' {Ar[^;C7'([/T ,(/+l)T ])] +J R z }. 
The association of the two cases gives 

Af[y M > N ;C?([0, (I + 1)T ])] < c^ 7-7 ' jl + Af[y M > N ; C^'([0, (/ + 1)T ])]| + ^. 



Since, for any t e [0,(Z + 1)T ], ||j/f' JV || By < |M|e y + M[y M > N ; Cj' ([0, (I + 1)T ])], one 
deduces 

A% M >VT 7 '([0, (/+1)T ])] < 11^11^+2^+247^ {l+JV[y M ' N ;C°>i'([0, (I + 1)T ])]} . 

To complete the proof, it now suffices to pick T such that 2c 1 x Tq~ 1 =1/2 and to set 

^ +1 = 2||V|| Sy +4^ + 1. 

□ 

3.3. Space discretization. This is the final step, that will lead us from y to y M,N . 
As in the previous subsection, we extend y M ' N on [0, 1] via (I38|) and use the notations 
r M,N ^yM,N,i introduced in Lemma [3.91 

Lemma 3.11. There exists a function C : (M + ) 2 — > IR + bounded on bounded sets such 
that if t p < s < tp+i < . . . < t q < t < t q+ i, with \t — s\ > 1/M, then 

11 M,m\\ < if i7' 

WVts ||By S M7+7 '-l \ l S \ ■ 

Proof. Thanks to the uniform control given by Proposition 13.101 one has 

II M,M,h\ 

\\Vts lls y 

q-2 

- II M,M\\ . 11 M,M\\ , || M.A/ || . \^ I , , 1-7' n M,M \\ 
S ||fly + lF stp Ha,, + lFt,t 9 _ 1 ||B y + cy 2^ F - tk+i\ Vt k+1 t k \\B 



1 l / J_ _ r y 

, \t - S| V 

° x ^\ M7-7' + M^+7'-i ( - Cl '*M^'-i 
where, for the last inequality, we have used the fact that 1/4<7'<1/2<7<1. □ 

Lemma 3.12. There exists a function C : (M + ) 2 — > M. + bounded on bounded sets such 
that if t p < s < t p+ i < . . . < t q < t < t q+ i, with \t — s\> 1/M, one has 

b tu dx u {F M -ld)fi{y u J|| By < M 2(j-y) I* ~ s l ■ 



k=p 

9-1 
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Proof. As Pm commutes with the semigroup, one can of course write 
fs tu dxtf (P M - Id) 

J s 

= X*f> M {P M - + (Pm - ld)A ts (X x ' l ' M 5f t (y M ' M )). 

From this expression, the uniform control given by Proposition 13.101 easily yields 
\\X^ M (P M - Id) f t (y?> M )|| Sy < c x \t - s\^' \\(P M - im{vf' M )h 



\t-s 



7 



while 



(P M -ld)A ts (X^ M 6f t (y M ' M ))\\ By 



< T7^II^(^' J ' M ^(/ /,M ))I| 67 < c x / 



□ 



We are now in position to prove the main result of this subsection, which, associated 
to Corollary 13.81 completes the proof of Theorem 12.41 

Proposition 3.13. There exists a function C : (IR + ) 2 — > M + bounded on bounded sets 
such that for any M G N*, 

sup \\y M -y M ' M \\ B y <C(||x|| 7 , W By )(||V-^/* y + T7^) ■ ( 46 ) 

ke{0,...,M} 7 I M7+7 j 

Proof. As in the previous subsection, we use the shortcut 

N[y M - y MM , CT' V (/)] := N[y M - y M > M ; C?(J, By)} + N[y M - y M ' M ; Cj'(I, By)}. 

Local result. Consider first an interval Jo = [0, To], with To a time to be precised at the 
end of this first step, and let s, t G [0, T ]. 

1 st case: if t p < s < t < t p+1 , then 

s(y M - v MM )ts = (¥% - x^ M p Mfl ( y ^ M ), 

hence 



\\S(y M - y MM )tsK < c 4 ,, x \t - s r < 



\t- s\ 



2 nd case: if t p < s < t p+ i <t< t p+2 , we go back to the previous case by noticing that 

WXl^-M „,Af,M\ II ^ 1 1 Jz-r-AJ „ M,M\ || i HXfrrrM „ M,M\ || 

\\o{y -y )ts\\B, < \\o(y -y K +1 |k, + \\6{y -y K+ lS |k,- 
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3 rd case: t p < s < t p+1 < . . . < t q < t < t q+1 with \t - s\ > 1/M. Then 

S(y M - y MM h 

's tu dx^ [f^)-P M n^ M )\+y^ m 

= fs tu d^ [Mrf?) - Mv?' M )] 

J s 

+ jf s tu dx^ (id -P M )fi(y^ M ) + y^ m - 

According to the two previous lemmas, one can assert that 

rt \t — si 7 ' 

: ;M (JA _ P. A f.(,,MM\ , . M,M,J r 



S tu dx\ M (Id -P M ) fi ( y M>M) + y^\\n, < C+, 



Besides, it is not hard to see that 



ii jT s tu dx^ im 1 ) - uy^ 1 )} k < 4 >x \t - «r' T^mr - y M > M ; c o ' v ao, m 

for some constant c]p x that we fix for the rest of the proof. 
Summing up the three cases, we get 



m M _ y u,M. C V ([0) Tq] . By)] < + c i ^'Mit 1 - y M ' M ; er' ([0, T ])]. 

In order to estimate J\f\y M — y M < M ; C°([0, T ], By)], it now suffices to observe that — 
y M,M = _ y M,M) g0 + Ss0 (ip - P M ^) 7 and so 

M[y M -y MM -A n \[^T })} 

< H-0 — ^^-011^ -K ^^^r 2 o^^T^-^'^V-[^ — 2/^'^; <5?^'([0, ^ ])]- 
Thus, pick T such that 2c^ x Tq~ j = 1/2 to obtain 

Af[y M -y M ' M ;C^'([0,T Q ])} < 2||^ - P M ^, + J^fc. (47) 

Extending the result: By following the same steps as in the local reasoning, we clearly 
get, for any rj > 0, 



Af[y M -y M ' M ;Ct([T ,T + V ];By)] 

< + cW^'^ 1 ~ y MM ; Ci°' 7 '([0, T + rj\, By)], 



which, together with (HT|) . leads to 
Af[y M -y M ' M ;Cf([0,T o + r,];B Y )} 

<2U- P M ^\\ By + + c\^'M\y M - y M > M ; ^'([0, T + v ], B Y )}, 
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and then 

T— Af „,M,M./>°a'l 



M[y M -y M ' M ;Cy([0,T Q + r,})} 



10 c 2 ^ x 



By taking r\ = T , we deduce 

Aay M -y M 'Vi°' 7 >,2T ])]<10 



20 cj x 



We repeat the procedure until the whole interval [0, 1] is covered. 

4. Rough case 



□ 



We now turn to the proof of Theorem 12.51 Thus, let us fix 7 G (g, 5), 7' G (1 — 7, 27], 
if) G By, and suppose that Assumptions (X2) 7 and (Fl) 3 are satisfied. We will follow 
(almost) the same steps as in the previous section: we first use pre-existing continuity 
results to reduce the problem to the study of the Wong-Zakai approximation y M , and 
then lean on a uniform bound for y M ' N to control the transition from y M to y M,N . 

Before we trigger the procedure, let us remind the reader with a few considerations taken 
from [S] on how to give sense to the equation under Assumption (X2) 7 . As in the Young 
case, the interpretation is based on an expansion of the regular equation: observe that 
if x is a piecewise differentiable process, then 

* s tu d&i f t {y u ) = xS7i(v.) + Xu'HfM ■ fi(v.)) + Jl 



'ts 1 



(48) 



where the operator-valued processes X x,t ,X xx ' l i have been defined by 
^S^dxiMis, with 



and Jf s :- 



Ml 



dr [f'iiVs + r(Sy) 

1 



+ 



^usVs 



fi{Vs)\ ■ (Sy) us 
S uv dx l v 5(fi(y)) vs 



cLuv dx^ v fj {ys) 



■fiiVs). (49) 



On top of the result of Lemma 13. 1\ one can here rely on the following extensions (we 
also anticipate on the sequel by introducing the additional process X ax ): 



Lemma 4.1 ([9]). The sequence of operator-valued processes 



X, 



ts 



Qitu dx u 



resp. X. 



X X ,JJ 

ts : = 



S tu dx^(Sx 



converges to an element X ax ' 1 (resp. X xx,lJ ) with respect to the topology of 

C2 +K (£(B a+K ,B a )) {a >0, kg [0,1)), 



Moreover, 



resp. C z 2 ~ ( - K (£(B ai B a+K )) (a G R, k G [0, 2 7 )). 

Af[X xx ^;C 2 2 ^ K (C(B a ,B a+K ))} < c a , K ||x|| 7 , 
M \XxM*",a -x xx ^-C^- K (£(B a ,B a+K ))} < c a , K v M , 
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and the same controls hold for X ax ' 1 in C^ +K (C(B a+K ,B a )). Finally, X ax ^ and X xx ^ 
commute with the projection Pn and satisfy the following algebraic relations: 

(SX xx ^) tus = X^(Sx% s , Xp* = - (Sx% s , (50) 

where X x ' 1 is the process given by Lemma \3.1\ 

Now, from a heuristic point of view, if we go back to the 7- Holder process x in (1481) . the 
expression (jl91) allows to identify J v as a B- valued process of order \i : = inf (37, 7 + 7') > 
1. This (partially) accounts for the definition: 

Definition 4.2. Let k G (0, 1) and i/j G B K . A process y : [0, 1] — > B K is said to be a 
rough solution of |?J) in B K if there exists two parameters \i > l,e > such that 

y = ^ and 5y - X^My) - X^{fl(v) " fM) e <2(B) nC5(B«). 
In accordance with the above considerations, one has in particular: 

Proposition 4.3 ([8]). If x is a piecewise differentiate process (resp. a standard Brown- 
ian motion) and if the initial condition ip belongs to B^ with r\ G (0, 1) (resp. r) G (|, 1) ), 
then the classical (resp. Stratonovich) solution of (T7|) is also a rough solution in B v . 

Remark 4.4. Let us go back here to the Young setting, ie when 7 > 1/2. In order 
to connect the above interpretation of (|7j) with the notion of solution derived from 
Proposition 13.4} observe the following equivalence: under the assumptions of Theorem 
13.5} a process y G CJ (By) is solution of ([7]) (in the sense of Proposition 13.41) if and 
only if y — and there exists ji > 1, e > such that 8y — X x > l fi(y) G C% (B) fl C|(By). 
Indeed, if y is the solution given by Theorem 13. 51 then, owing to the contraction property 
([20]), 5y - X'tfiiy) = A(X x ' i 5f i (y)) G C 2 7+V (B) n C^(B y ). On the other hand, if 5y - 
X x ^f t (y) e C$(B) n CI(B Y ) and z is defined by z = ^5z = X^f^y) + A.(X x ' i 6f i (y)), 
one has S(y — z) G C^B), with jl = inf(/i, 7 + 7') > 1. As yo = Zo, this easily entails 
y = z. 

4.1. Previous results. With the above definition in mind, the main result of [S] can 
be summed up in the following way: 

Theorem 4.5 Q£j). Under the assumptions of Theorem \2.5[ Equation admits a 
unique rough solution in By in the sense of Definition \4 . 2\ Moreover, if y (resp. y) is 
the rough solution in By of ([?]) associated to a process x (resp. x) that satisfies (X2) 1 , 
with initial condition ip (resp. ip) in By, then 

M[y - y; (*([0, 1]; By)} < C (||x|| 7 , ||x|| 7 , ||^|| By , ||^|| y ) { ||^ - ^|| By + ||x - x|| 7 } , 

(51) 

for some function C : (M + ) 4 — > M + bounded on bounded sets. 

As in the Young case, denote y the Wong-Zakai solution of (EJ), which corresponds 
to the classical (or equivalently rough) solution of the equation when x is replaced with 
x 2 . The continuity result (15T|) allows to control the transition from y to y M : 

Corollary 4.6. Under the assumptions of Theorem \2.5[ there exists a function C : 
(M + ) 2 — > R + bounded on bounded sets such that for any M , 

sup \\y th -yf k hy < C(\\x\\ y , \\ip\\ B ,)v 2 m. 

k£{0,...2 M } 



IS 
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Now it is worth noticing that the time-discretization of the equation has been analyzed 
in [8], too. In other words, we already know how to control the difference between y M 
and the process y M generated by the intermediate Milstein scheme: y$ = ip and 

iC = s tk+ltkV ? k + x^&Myg) + xij k « (/;«) ■ /,«)) , (52) 

where t k = tf = To express this result, let us denote (n M ) the sequence of dyadic 
partitions of [0, 1], and introduce the two processes 

K%:=(6y M ) ts -X? s -7 4 (y^) , JM:=KX-XZ x « (/^f ) • £(</f)) , 

for every s < i 6 I1 A/ . For sake of clarity, we will also appeal, in the sequel, to the 
discrete versions of the generalized Holder norms introduced in Subsection 12.51 Thus, 
for any M e N, we denote [a, 6] M := [a, b] H U M and 

M[h- c*at;, t»}„, B a , p )} ■.= sup IIW^II^ , 

tp<s<t<t™ \t — S\ 

s,teu M 

We define the quantities 

A^C^MIm;^)] , Ar[.;C 2 A ([a,6] M ;^, p )] , Af[.; C 3 A ([a, 6] M ; B a , p )], 
along the same line. 

Proposition 4.7 ([H])- Under the assumptions of Theorem \4-5[ for every 

< p < inf ^7 + y' - 1, 7 - 7' + i 

t/iere exists a function C = Cp : (IR + ) 2 — >■ R + bounded on bounded sets such that for any 
M, 

||-M A/ 1| ^ C(||x|| 7 , NjKj 

sup ||y tfc - y ifc || B < (53) 

fc=0,...,2 M / 

where y is the process generated by the intermediate Milstein scheme Moreover, 
there exists another function C : (M + ) 2 — > M + bounded on bounded sets such that the 
following uniform control holds: For every M , 

Af^-C^lOAjM-By)} +A% A/ ;C7([0,1] M ;£)] + Af[K M ; C 2 2 ^([0, 1] M ; B)\ < c x4 ,. (54) 

where c x ^ := C"(||x|| 7 , 



It now remains to study the transition from y M to y M,N , which is the purpose of the 
two following subsections. 

4.2. A uniform control. The aim here is to exhibit a uniform estimate for y M ' N , to 
which we will extensively appeal in the next subsection. As in the time-discretization 
procedure, the two following processes will play a prominent role in our reasoning: for 
every M, N and every s < t G I1 M , define 

R M,N ;= {h M,N )u _ xf'PxMy^), 

J^ N := K»>» - xf* 2M "P N (f>(yf> N ) ■ P N fAvf' N )) ■ 
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Proposition 4.8. There exists a function C : (M + ) 2 — > IR + bounded on bounded sets 
such that for every M, N, 

Af[y M ' N - d?([0, B)] + Af[y M,N ; C°([0, 1] M ; By)] + Af[K M > N ; C 2 2 ^([0, 1] M ; B)] 

<C7(||x|| 7 ,||^||^). (55) 

Proposition 14.81 is actually a spin-off of the following successive controls on J M,N : 
Proposition 4.9. Fix e, fi such that 

7 + 7'>yU>l , 7- (V - > £ > °' 

There exists two integers M = M (||x|| 7 ), iVo = JVo(||x|| 7 , HV'IIb /) onc ^ time T = 
T (||x|| 7 ) > 0, T G n Mo ; snc/i ffcat /or every M > M , N > N , for any k, 

A/y^; ^([fcT , (* + 1)T A 1] M ; B)] < 1 + \\y^\\ B y (56) 

and 

Af[J M ' N ; CI(pT 0) (A + 1)T A 1] M ; ft/)] < 1 + II^TlUy (57) 

The proof of Proposition 14.91 resorts to the following technical lemmas: 

Lemma 4.10 Q8j). Let e > and \i > 1. There exists a constant c = c EtfJi such that for 
any M and any process A: S 2 By satisfying A t M ^.m = for every k G {0, . . . , M— 1}, 

\\A ts \\ B <c\t- s\»M[6A; C£([s, tj M ; B)} (58) 

and 

||A s || By < c{|t - S | £ + I* - *r 7 '} {^[M;C^([ S ,t] M ;B)] +^[M;C 3 e ([ S ,t] A/ ;B y )]} . 

(59) 

/or any s < t G IT . 

Lemma 4.11. There exists a function C : 1R + — > M + bounded on bounded sets such that 
for any M and every s < t G IT M , 

Af[6J M > N ;C»(ls,q M ;B)} < c x \t - s\^'^ 

J Af[y M,N ; C?([a, By] + Af[y M ' N ; tj M ); B] + Af[K M ' N ; C 2 \{s, tj M ); B] 

\ + JV 1 - 2 ^^; *]*,; B y )] 2 ) (60) 



and 

Af[6J M > N ; Q(ls, t\ M - By)} <c x \t- s \^'-^- e 

{l+N[y M ' N -A{{sAM-,By)\} {l + N 1 ^'Af[y^ N ;C 1 (ls,tj M ;By)] 2 } , (61) 

where c x := C(||x|| 7 ). 
Proof. See Appendix. 

□ 
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Proof of Proposition \4-9[ For sake of clarity, we write here x for x 2M . 

Step 1: k = 0. This is an iteration procedure over the points of the partition. Assume 
that both inequalities ([55]) and (J57|) hold true on [0, t^J M and t^ +l < T (T will actually 
be precised in the course of the reasoning). Then, for every t G [O^^Jm, one has 



WVt ■ II 



M,N 

< \\Jw' N K + H^I| By + wxgPyfiWWz, + \\x^p N (m ■ p N fM))\K 

< i + 2|mi b . + Cx {t^'-^um\B 1/2 + t^'-^wm ■ p N mh,' 



where 77 G (0, |) is picked such that 27 > 7' — 77. Now ||/i(^)||B 1/2 < c|l + ||^||b 7 ,| an d 
as in (JHZD and 

r 1 f UWb 

\\m-PNfm^<c\i + \m B j\i + 



AT27'" 



hence 

r i f IIV'IIb 



VJ I ]\T2 7 '-1 
1 1 II^Hs , 

At this point, let us introduce an integer N ' such that ^l.i^y-i < 1; so that for any 
N > Ni*, M[y M ^-Am^ hi; B Y )} < 2c\ {l + ||^|| By }. Besides, if a < t G [0,tf ] M , 



II (Sy M ' N 



ts\\B 

< Uf N h + K/PnMv? > N )h + U^PMiv^) ■ P N fM UN ))\\B 

< c,|t-s| 7 {l + ||^|| By }, 

and 

\\K t f N h < Uf N \\B + Kr j p N (fl(y^ N ) ■ Pn/Av^Mb 

< C |t- s r{i+iHi By }. 

By using ( |60|) . we deduce, for any N > Nq 1 (remember that t^ +l < T ), 
Ar[5J A ^;Cr([0,tJ +1 ] M ;i3)] 



AT27'- 



clUP 



< cjr'"'1i T yi, Ui 



and in the same way, according to (I6T1) . 

Ar[5J^;C3 £ ([0,tJ +1 ], /; ^)] < c,Tr (V -^- e {l + |H| By } J 1 - ^"^"^ 



iV 2 7'-l 
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Let us fix an integer N Q ' > N ' such that N2l i-i < 1- Then, thanks to ( 158J) and (J59 



N 2 r 

(applied to A = J M ' N ), we obtain, for any N > N ' , 

U[J M > N -Amtf + M, B)\ < 4 + rr (7 '" l)_E } {i + W» v }. 

and 

^[^^([O^ilMjfi/)] < cl {T^ + Tr^- 6 } {l + ||» y } . 

Consider now a real T ' > such that c 3 x {(7^+^ + (T^)^^'-^- £ | < 1 and let 

M = M (||x|| 7 ) an integer such that 1/(2 M °) < Tq. We fix T in the non empty set 
(0, Tg)nll Mo so as to retrieve the expected controls, namely: for every M > M ,N > Nq 



N[J M ' N -Amtf + iU-M < 1 + IHly, ^[J M ' N ;CI(l0X + xiM;By)} < 1 + , 

which completes Step 1, that is to say the proof of (15E|) and (1571) on [T , 2T J M . 

Step 2: k — 1. We henceforth fix M > M . With the same arguments as in Step 1, we 
first deduce, if both controls (loT)j) and (lo"T|) are checked on [T ,tf] M (with tf +l < 2T ), 

f \\ V M > N \\ 2 1 

Ar[^; Cl ([T ,tf]M;^)] < 4 {l + h% N \K} 1 + ■ (62) 



Remember that for any N > Nq' 1 , WvtI^Wb^ — 2c* |l + H^ll-By }• Consequently, 



we 



2 1 12 (24{l+||^||s ,}) 2 

introduce an integer N ' > N ' such that — ^a.ipy-i — < 1; an d dS2J) entails, for any 

N > Nq' 1 , 

Af[y M ' N ;mT ^j M ;B y )} < 2c\ {l + H^H^}. 
Then, with the same estimates as in Step 1, we get, for any N > Nq' 1 , 

1 ( c 2 \\v M,N \\ 2 1 

W^mn^jM-M <Ui+ ''^JM {i + n^iia,} 

and this control also holds for M[J M ' N ; Cf([T , tf +1 \ M ; By)}. Introduce finally iV 2 ' 2 > 

Nq 2 ' 1 such that Ci{2C ^^ B j' }? < 1 to get (jSE} and (ETJ) on [T ,^J M (for iV > iV 2 ' 2 ), 
which completes the proof of (155]) and (15T1) on [T , 2T ]m- 

We repeat the procedure until Step L, where L = L(||x|| 7 ) is the smallest integer such 
that LTq > 1. 

□ 

Once endowed with the estimates of Proposition 14.91 the proof of Proposition 14.81 
follows the same lines as the proof of Theorem 2.10 in [8]. For sake of conciseness, the 
reader is referred to the latter paper for further details on the procedure. 
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4.3. Space discretization. We are now in position to compare y M with y M,N . To this 
end, let us introduce the intermediate quantity: for every s < t G I1 M , 



+ M[y M - y M > N ; CJds, tj M ; B)] + M[K M - K M ' N ; C 2 2 ^([ S , tj M ; B)]. 



Lemma 4.12. For any A G (0,7 + 7' — 1), there exists a function C — C\ : (1 
bounded on bounded sets such that for any M, N, for every s<t6 I1 M ; 



AT[5(J m -J m ^);C 3 ^- a ([ S ,^m; 8)1 < <| |t - s| A jV[y M - y M ^; Q(M]m)] + 



and 



N 2\ 



1 1 



M[8{J M -J M > N )-Cr\{sAM-,B,,)] < ||t - s| A A% M - y A ^; Q(M]a/)] + ^ j 
w/iere := C(||x|| 7 , ||^|| By ). 

Proof. See Appendix. □ 

Proof of Theorem \2.b\ For sake of clarity, we write here x for x 2 *' . Consider a time 
Ti G [0, 1] M - For any t G [0,Ti] A /, one has 



M M.N 

Vt -Vt 



Thanks to Lemma T4. 101 (applied to A = J M — J M,N ) and Lemma T4. 121 we already know 
that 

\\J% - Jto' N K < {TiN[y M - y M ' N ; Qmm M )} + n~ 2X } . 

Then 



< c 



p 



N- 



1 

-2A 



and with similar arguments 



< c x ^\\\^-P N iP\\ B .+N- 



-2A 



so that 

M[y M -y^-AmmM-B,,)] 

< {T>M[y M - y M ' N ; Q^O^m)] + U~ P N ^, + N~ 2X ] 
Let us now analyze (in B) the decomposition: for every s < t G [0,Ti] A /, 

% A/ - y M ' N )ts = [xtfhiyf) - x^p N f l ( y ^ N )} 



[x x t r (myf) ■ fM)) - x?r j p N (//^ . p N fM> N ))] 



J, 



M 
ts 



J, 



M,N 
ts 
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According to Lemmas 14.101 and 14.121 

H4f " Ju' N h < \t - s\* {T7'-W[^ - y M > N ; QQO^Jm)] + AT 2 *} . (63) 
Moreover, 

\\X^h{yf)-XtfP N f i {y^ N )h 



< c x \t-sr \\yf - yf> N \\ B + ||Xf/(Id -P N )Myf' N ) \\ B 



N -2\ 



< c x4 , \t - s\" {^M[y M - y M > N ; QdO,^)] + U ~ Pi 
and this kind of argument leads to 
M[y M -y^-AmUM-B,,)} 

< {nM[y M - y M > N ; Q([0, T^m)] + U~ Pn^\\b v + N~ 2X 

Finally, 



K ts ~ K ts 



T M _ jM,N 
J ts J ts 



[xtr (fM 1 ) ■ ay" 1 )) - x;:-'i\ (m y ?> N ) ■ jw*(i^))] + 

and thanks to (1631) . this decomposition easily allows to conclude that 

N[y M — y M ' N i Q([o, Ti]j\/)] 

< 4,* {r/-W[^ - y M ' N ; Q([0, TJm)] + U - P^llsy + A~ 2A } - 

Let T* > and M G N such that 4^(T*)^' = \ and (0,T*) n n M ° ^ 0. The time 7\ 
is now fixed in (0, Tj*) D n M ° so as to retrieve, for any M > M , 

M[y M - y M > N ; Q([0, ft]*)] < 2c^{u ~ Pjv* y + N~ 2X ] . 

It is readily checked that the same reasoning (with the same constants) holds on any 
interval \kT\, {k + l)Xi A 1}m and leads to 

M[y M - y M ' N ; Q{{kT u (k+)T 1 A 1} U )] < 2c^ { \\yjf Tl - yjgf || By + A" 2A 
As ft only depends on x and it follows from a standard patching argument that 

Af[y M - y M ' N ; ^([0, &/)] < c,, v ,iV" 2A , 

which, together with the results of Corollary 14.61 and Proposition 14 .7} completes the 
proof of ffTB^) . 

□ 



5. Appendix A 

Let us go back here to the technical proofs that have been left in abeyance in Section 
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Proof of Lemma ^777 For sake of clarity, we write here x for x 2 and y for y M,N . One 
has 

(5J M ^) tus = X^P N 5(f\(y)) us - X^(8x% s P N UKVs) ■ P N fAVs)) 

+ X?:> ij P N 6(fl( y )-P N f j (y)) us , (64) 

which easily entails 

(SJ M ' N ) tus = 1%? + II&T + + IV t T, (65) 

with 

Cs N ■= XgPN (f dr f>(y s + r(Sy) us ) ■ K^ N ) , 

II t t N := XgP N (J dr f'{y s + r(5y) us ) ■ [a us y s + X^Pn^Vs)] \ , 

III^ N := XgP N (j dr [ti(y s + r(5y) us ) - //(„.)] (5x\ s ■ Psffa) 
IV t ^ N := X^P N S (My) ■ P N fj(y)) us • 
First, \\l£F\\ B < c, \t - «r II^IIb, and 

Unfile < c x \t-u\^\\u- S r' hsh^ + lu-s^WMy^ 



3y ' I"- I |U«\ws/||Oi/2 

< c.|t-sr +y {i + w By }. 

Then 

II^^IIb < c a .|t-s| 2T ||(«Jy) 1ia ||B||Pjv/i(y,)||i,<» 

< c x \t- s\ 2 ^ {\\(Sy) us \\ B + \u- S | 7 ' \\y s \\ By } {1 + ||(Pv - Id)/^)^} 

and 



N 2 y 

c 

< 



py —{i + \\ys\\^\- (66) 

Finally, 

\\IV t Th < c x \t-u\^\\(5y) us \\ B {l + \\P N f 3 (y s )\\ L ~} 



\\y II 2 

< c x \t-u\^\\{8y) us \\ B ll ' 8 B -<' 



Going back to (1551) . those estimates yield (15U|) . To get (15T1) . we resort to the decompo- 
sition ( I64p and observe that (for instance) 



1/2 



< Cs | t _ sr (7'4)|i + MBy | 
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and for any 77 G (7' — 7, |), 

\\X^(8xi) us P N (ti(y s ) ■ P N f j (y.))\\ Bl , 
< c m \t-ur«-4\u-8F\\fl(y.)-P N f j (y. 

| 2 7 -( 7 '- 7? ) 



{l + || ? / s || By }{l + ||P J v/ i ( 2 / s )||L0 O }, (67) 



where, to get the last estimate, we have used the property fl2T|) . Together with fl66|) . this 
leads to ( ISTj) . 

□ 

Proof of Lemma \4.12\ Observe first that SJ M can be decomposed as in (IMj) or as in ([65j). 
by suppressing in both expressions the projection operator Pjy. In order to estimate 
\\5(J M -J M > N ) tus\\i3, we re ly on the decomposition fl25]) and its equivalent for J , with 
J M instead of I M > N ^ etc. Write for instance 

jM _ jM,N 
tus tus 

= X% Qf * dr [X(y? + r(5y M ) us ) - f'^ N + r(6y M > N ) us )] • K% 
( ^ dr fl(y^ N + r(5y M ' N ) us ) ■ [K% - K% N ] 



+X-(Id -P N ) (J^ dr + r(6y M > N ) us ) • =: + J« + /£>. 

Owing to the uniform estimate (151| and the continuous inclusion By C one has 
first 

< c x>iJ \t-s\^M[y M -y M > N ;Q(ls,tU)}. 
Then clearly ||/ t ( 2|| B < c x \t - s\^ M[K M - K M < N ; C?{{s, t\ M] B)] and 

\\I%sh < c x \t-ur X \\(ld-P N ) (J\r f'(yr> N + r(5y M > N ) us ) • K^\\ B _ X 

< c x \t-ur x N- 2X \\K% N \\ B < Cx ^\t-s\^-\ 

where, for the last estimate, we have used the uniform control (|55|) . The other terms 
7, II, III of (1651) can be handled with similar arguments. Let us only elaborate on 
the estimate of \\X™' ij P N (fl(y^' N ) ■ (Id -P N )8fj(y M ' N ) us ) \\ B , which may be a little bit 
more tricky. Indeed, one must here appeal to the property fl22|) to get 

WPn (mVu' N ) ■ (Id-P^Sfjiy^^Ws 

< c x \t - u\^- x wm y ™> N ) ■ (H -pnW^uib^ 

< Cx \t- u \^~ x \\my^ N )\\ B ,, im-p^sf^ui^ 

< c St 4\t-u\* r - x N- ax \\6f j (y M ' N ) u .\\B < c x ^\t-s\ 3 ^ x N~ 2X . 
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As far as \\5( J M - J M > N ) tus\\B i is concerned, one can start from the decomposition (1641) 
and observe for instance that 

< [My?) - h{y^ N )} Ha, + ||X^(Id -P N )Uy^ N )\\ Bl , 

< c x \t - «r II f dr fl(y^ N + r(yf - y^)) ■ (y™ - y^ N ) || Sy 

Jo 

+c x \t-ur x \\(ld-P N )f, t (y^ N )h lf _ x 

< c^{\t - s r ibf - y ^ N h^ + it - s r x n-^} . 

The other terms steming from can be estimated along the same lines. □ 

6. Appendix B: Implementation 

We would like to conclude by insisting on the simplicity of the two algorithms (TTOT) and 
(TTTj) as far as implementation is concerned. To this end, we focus on the case x = X is a 
fBm with Hurst index H £ (1/3, 1) and x is its linear interpolation. As pointed out in 
Section [21 we know that x satisfies Assumption (X2) 7 (and accordingly Assumption 
(Xl) 7 )for any 7 £(§,#). 

6.1. Young case (H > 1/2). The objective here is to implement the Euler scheme 
(TTOj) . Remember that we have fixed a basis (e n ) of L 2 (0, 1) made of eigenvectors of A. 

By setting Y t ^' M ' 1 = (y/^ ,M , e^y, one has, for any I £ {1, . . . , M}, 

A/f m 

y m,m,i = e - k/ M Y M,M,i e . Xl/M] J2(5X*) tk+ltk (h{Y t f M ), e) . (68) 

1 i=l 

The following Matlab code is a possible implementation of this iterative procedure, 
for which we have taken m = 1 and 

^(£) = isin(7r£) + |sin(37r0 (£e[0,l]), fk { x ) = tSLzA {x G M ). (69 ) 

The parameter /c is meant to vary so as to observe the influence of the perturbation. 

The procedure more precisely simulate the evolution in time of the functional-valued 

process Y M,M . At each step, the Fourier coefficients (^fi(Y t A k 4,M ),ei^ are computed by 

means of the discrete sinus transform function dst (and its inverse idst), according to 
the approximation formula 

J ° n=0 

As for the fBm increments, they are computed via (an approprietly rescaled version of) 
the Matlab-function wfbm, which lean on the decomposition of the process in a wavelet 
basis, following the method proposed by Abry and Sellan in [TJ. Let us finally point out 
that the action of the semigroup is likely to be qualified by turning the heat semigroup 
S A into St = S^j., for some parameter k. The theoretical study contained in Section [3] 
remains valid for the modified system, of course. 
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function [ l] = eigval (N) 

1=[] ; for i=l:N, 1 ( i ) = ( pi * i )" 2 ; end 

function [S] = semigr (M, N, 1 , kappa) 

S = [];for i=l:N, v ( i )=exp(- 1 ( i ) " 2/ ( kappa*M) ) ; end 

function=simulyoung (H,M,N, k , kappa) 
l=eigval (N) ; S=semigr (M,N, 1 , kappa ) ; 
X= ( 1 /M) * H* wfbm (H ,M+ 1 ) ; 
A= [1/2, 0,3/5, zeros (1 ,N-3)]; 
u=zeros ( 1 ,N) ; fy=zeros ( 1 ,N) ; 
for i=l:M 

u=dst(A(i ,:)); fy =0.5* idst (k*(l -u) . / ( 1 +u . " 2 ) ) ; 
A(i+l,:) = S.*A(i ,:) 

+ ((kappa./l).*(l-S))*M*(X( i+l)-X( i )) .* fy ; 

end 

E=[];for j=l:M+l, E(j ,:) = dst(A(j , : ) ) ; end 
plot (linspace (0 ,1 ,N+2) , 

[0 , dst ([1/2 ,0 ,3/5 , zeros (1 ,N-3)]) ,0]); 
F(l) = getframe ; for p = l:M 

plot (linspace (0 ,1 ,N+2) ,[0 ,E(p + l ,:) ,0]) ; 

hold off; 

F(p+l) = getframe ; end 
movie (F, 1 , 2 ) 



Figure IBTTI corresponds to simulations of the process t h- > Y t ' (|) for different values 
of the parameter k (k = 1,5,20,50), with the same realization of a fBm with Hurst 
index H = 0.6. The above-described parameter k has been taken equal to k = 100. 



6.2. Rough case (H G (1/3, 1/2]). By projecting Y ' onto e/, one retrieves here the 
iterative procedure: 



= e-^"Y t f* + 2 -{l- e-^"}j2(SX%^ (/,«*»),«) 

1 i=l 
m / ftk+1 \ 

+ (2 M f (^% k+ltk (5X% +ltk / e-^^du (u - t k ) 

(PNf^ N )-f'(yf> N ),e^. (70) 
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M k=5 
-0.2 1 1 1 1 1 1 1 1 1 1 1 -0.2 1 1 1 1 1 r- 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 ' 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




Figure 1. Influence of the perturbation term through the observation 



of the path t H- Y t RI,N (~), for different values of the parameter k in ([69 
(k = 1, 5, 20, 50). Here, M = N = 1000, H = 0.6. 



The computation of the Fourier coefficients (fi(Y% ,N )-> eij can be implemented with the 
discrete sinus transform, as in the Young case. As for the computation of 

it can be achieved with the same idea, starting from the approximation: 

( M,N\ nil M,N\ 

p Nfj{y tk ) ■ fi(vt h )> e 

^ N N N 

N2 E E E e < (a?) ^ (a?) e ? Q fi (y" ,N (a?)) u (vi 



M.N 



n=0 p=0 m=0 

Those considerations easily lead to the construction of an algorithm for flTOj) . 
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